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Abstract 

We use here a particle system to prove a convergence result as well as a deviation inequality 
for solutions of granular media equation when the confinement potential and the interaction 
potential are no more uniformly convex. Proof is straightforward, simplifying deeply proofs 
of Carrillo-McCann-Villani [CMCV03, CMCV06] and completing results of Malrieu [Mal03] in 
the uniformly convex case. It relies on an uniform propagation of chaos property and a direct 
control in Wasserstein distance of solutions starting with different initial measures. The deviation 
inequality is obtained via a Ti transportation cost inequality replacing the logarithmic Sobolev 
inequality which is no more clearly dimension free. 
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1 Introduction 

Our main goal will be to deal in a probabilistic way with the following nonlinear equation 

— = V • [Vn + uW + uVW * u] (1) 
at 

where u{t, x) is a time dependent probability measure, * denotes the standard convolution operator 
and V and W are two convex (at infinity) potentials. 

This equation arises (in dimension 1) in the modeling of granular media as follows: consider many 
infinitesimal particles colliding inelastically. With a correct renormalization between the frequency 
and the inelasticity of the collisions, u{t, x) turns out to be the velocity of a representative particle 
(among an infinity). The potential V represents the friction and W the inelastic collisions between 
particles with different velocities. Note that the particular case ^ = and W{x) = is of 
special interest and we refer to Benedetto-Caglioti-Pulvirenti [BCP97] for the physical issues (also 
see [BCCP98]). 

Once the problem of existence and uniqueness is tackled, one major problem in this equation is 
the behavior at infinity: existence of a stationary measure and speed of convergence towards this 
stationary measure or even distance between two solutions starting at different points. It has been 
studied in parallel by Carrillo-McCann-Villani [CMCV03, CMCV06] and Malrieu [MalOl, MalOS], 
under various assumptions on the potentials V and W , using analytical and probabilistic approaches 
respectively (also see [BRV98, BRTV98] for one dimensional particles). We will consider here the 
probabilistic approach and will recover and generalize slightly results of [CMCV03, CMCV06] as 



well as give a quantified probabilistic approximation of the stationary measure of the granular media 
equation. It is worthwhile noticing that the analytic methods in [CMCV03, CMCV06] cover a much 
larger spectrum of non linear p.d.e's (like the porous medium equation), for which the probabilistic 
approach remains to be written. 

As physical interest (and in fact where the main mathematical difficulty resides) dictates the friction 
term to vanish, we will consider the following two sets of general assumptions: 

Definition 1.1 • We say that W satisfies the set of assumptions (A) if 

Al. the friction term V = 0; 

A2. W is symmetric, i.e. W{—x) = W{x); 

A3. VW is locally Lipschitz with polynomial growth, i.e. for some m 

Var, y G \VW{x) - VW{y)\ < C{\x - y| A 1)(1 + I^P + IvD- 

We also assume that the second derivative ofW has a (m) polynomial growth. 
A4. W is the sum of a compactly supported function and a uniformly convex function. 
This last assumption entails some uniform convexity at infinity property, namely that 
there exist positive C, X, such that 

Va;,y, {x - y).{VW{x) - VW{y)) > \\\x - yf - C . (2) 

• We say that V and W satisfy the set of assumptions (A') ifV is uniformly convex at infinity 
(i.e. satisfies (A4)), W is convex at infinity (but not necessarily uniformly) and symmetric, 
and both satisfy (A3). 

As we will sec later, these assumptions are not sufficient to get good properties for large time of the 
granular media equation and will be replaced by some "strict convexity except for a finite number 
of points" property. 

Before further discussing the result, let us first present the probabilistic approach of this problem. 
The probabilistic interpretation is to consider a Markov process {Xt)t>o, which law at time t is u. 
It is the solution of the nonlinear S.D.E. 

r dXt_ = V2dBt - VV{Xt)dt - VW * ut{Xt)dt 
\ C{Xt) = ut dx. 

We wish here first to give sufficient conditions ensuring both existence and ergodicity of the solution 
of the nonlinear S.D.E. and in a second time to provide a way to simulate this law at each time t 
with some Gaussian confidence intervals independent of time t. These two goals will be carried out 
through the extensive use of some (linear) particle approximations, i.e. (X^) solution of 

= V2dBl - VV{Xl''')dt - i Ef=i VTy(X^'^ - X/'^)di i = l,...,N 
X^'^ = X^ i = l,...,N. 

If we suppose that F = 0, we may assume that the center of mass is fixed, and without loss of 
generality, set to 0, i.e. EXq = 0, indeed it is easy to remark that Xt — ^Xq satisfies equation (3). 
This assumption obliges to introduce the projected system onto the set X^ja;^ = (see Section 2), 
and introduces some intricacies. 

Under the strong assumptions of (global) uniform convexity of one of the potentials (the other one 
being globally convex but non necessarily uniformly), Malrieu [MalOl, MalOS] successively fulfilled 
these two goals using extensively the so called Bakry-Emery criterion ensuring that a Logarithmic 
Sobolev inequality holds independently of A'^. He also provided in the same time asymptotic behavior 



(3) 



in time and concentration inequality, enabling him to recover part of the results of Carrillo-McCann- 
Villani [CMCV03] (see also Bolley-Guillin-Villani [BGV05] for a strengthened deviation inequality 

for almost quadratic potentials). 

The assumption of uniform convexity is however far too strong for applications (preventing for 
example to consider the case of cubic interaction potential). It can be removed (see [CMCV03]) for 
a one point degeneracy where the authors obtain various rates of convergence. Another approach 
is proposed by Carrillo-McCann-Villani in the subsequent [CMCV06] paper introducing tools of 
contractions in L^-Wasserstein distance for length space enabling them to ensure a convergence in 
Wasserstein distance of the solution of the granular media equation. 

We will see here that, with much simpler tools, we can recover their result and may provide also 
Gaussian confidence bounds for the approximations of the granular media equations. It is worthwhile 
noticing that, though the invariant measures of the particle system still satisfy some log-Sobolev 
inequality, it is difficult to obtain dimension-free estimates on the log-Sobolev constants (as for 
nearest neighbors interactions models in Statistical Mechanics). This prevents us to use this kind 
of approach. 

Our main tool in order to get the Gaussian concentration inequality we need, will thus be a trans- 
portation cost-information inequality: let say that G Ti(C) if for every probability measure v 



where Wi is the usual Wasserstein distance and Ent the Kullback information or relative entropy. 
In the final section we prove that, under one of the hypothesis (A) or (A'), the law of the particle 
system at time t satisfies a Ti{CN) inequality, for some C that does not depend on N nor on t. 

To complete the proof of the concentration inequality in the final section we need some "uruform 
(in time) propagation of chaos". It seems difficult to obtain such a result only assuming (A) or 
(A'). Indeed if the potentials are non attractive in some (bounded) region, the situation becomes 
the classical mean-field one where propagation of chaos is controlled on finite time intervals only. 
That is why we have to reinforce our assumptions. 

Our condition below is inspired by the work of Carrillo-McCann-Villani [CMCV06]: 
we say that condition C(A, a) holds if there exist A,a>0 such that for any < e < 1, 



Remark that this condition implies convexity {if \x — y\ < 1/2 choose e = |a; — y|, otherwise choose 
e = 1/4) but is weaker than uniform convexity (which is true with a = 0). Typical examples are 
polynomial potentials of the form W{x) = which satisfy C(A, a) for some positive A. In fact 

this condition is related to the rate at which strict convexity is lost at points. 

Condition C(A,a) allows us to prove the required "uniform (in time) propagation of chaos" and 
then to obtain the rate of convergence of ut to the limit Uoo in Wasserstein distance W2- These 
results are obtained in sections 3 and 4. 

To complete the description of the paper, the next section contains some useful estimates, and a 
complete proof of existence and uniqueness for all the equations we have considered ((1) and (3) in 
particular). The proof is quite natural (using the particle system). We did not find such a proof in 
the literature except in dimension 1, and as a matter of fact the set of hypotheses we need, slightly 
differs from the ones used in the analytical literature. 



2 The particle system and solutions of the nonlinear SDE 

In this section we study the particle system and the nonlinear SDE. The first discussion is about 
case (A) when V vanishes. 



(5) 



Vx, y e M'^, {x - y).{VW{x) - VW{y)) > Ae'^Qx - y\^ - e^). 



(6) 



Consider once again the particle system (X^) solution of 



(7) 



Here Xq are i.i.d. variables, with common law /^o- 
As noted by Malrieu, when F = 0, the direction {{v, ■, , v)} is quite singular for the particle system, 
and we then introduce its orthogonal hyperplane = |x G (M*^)^; X^^^ = o|, and consider the 
projected particle system 



1 ^ 



The process {Y^ )^^q thus verifies the system of SDE's 

/ dy/'^ = V2dBt - 4 Ef=i dBi - Ef=i vw^(y/'^ - i = i, n 



(8) 



hence is a diffusion on {M}. Note that E(y/'''^) = 0. 



We assume that W is satisfying the assumptions (A) in Definition 1.1. Let us first remark that 
this new particle system satisfies globally the "convexity at infinity" property. We have to verify, 
denoting b{x) = {b^{x), ...,b^ {x}), 6* with values in M'' and 



1 ^ 



that condition (2) holds true. Note that, with a the (constant) matrix diffusion coefficient of 
{Y^'^)l^(.i<J^ , we have < c for some positive c, and, as W is symmetric and x,y G M. 

1 ^ 

{x-y)- {b{x) - b{v)) = -— J2 i^i - - iVi - Vj)) ■ C^Wixi - Xj) - VWivi - %■)) 

1 ^ 

CN 



2N 



-X\\x-yf + 



2 ' 



where the latter equality follows as x,y G Ai. 

This remark allows to prove existence, uniqueness and non-explosion for the solution of the SDE 

(8). Indeed recall that a sufficient condition for all this to hold is the following: 

there exists some i/j such that il^{x) — +oo as \x\ — +oo and + b.Vt/j is bounded from above. 

'0 is a kind of Lyapounov function. Here we may choose ip{x) = |xp according to the previous 

remark applied with y = 0. In this case non-explosion holds as soon as /^od^^P) < +oo. 

The same argument can be used with (A') if W is everywhere convex. 

For the initial system (7) the only thing to show is non explosion. It will follow from Proposition 
2.1 below. 



2.1 Moment controls for the particle system 

In order to prove tightness of the empirical measure of the particle system, and for later use, we 
first prove some controls on moments. 

Proposition 2.1 IfW and V satisfy (A'), and //o(|a;p) < +00 the solution to (7) is non explosive. 
Furthermore there exists some K > such that for all i, 

supE|X;'^|2 < iiq{\x\^)+K. 

IfW satisfies (A) a similar result holds with Y^'^ in place of Xl'^ . 
Proof 

< Using Ito's formula up to the infimum between t and the exit time of a large ball (let T be this 
infimum) we have 

TV N N „T 

= + 2dNE{T) - 2E^ / Xl'^ .VV{Xi'^)ds 

1=1 i=l i=l 

-4^E / Xl'''.VW{Xl'^-Xi''^)ds 
i,j=l -^0 

N N T 

= E^|X^'^p + 2dArE(r)-2E^ / Xi'^.VV{Xi'^)ds 
1=1 i=i "'o 



We may use condition (2), go to the limit with respect to the radius of the ball (hence replace T by 
t) and obtain the finiteness of the quantity. This implies non explosion. 

Furthermore, denoting v{t) = E^^j^ differentiating (see the proof of the next Proposition) 

and using our conditions (in particular x.VW{x) ^ Av^la^P — C), we get 

v'{t) < -2Xvv{t) + 2dN + 3CN. 

Gronwall's lemma and exchangeability conclude the proof in case (A'). 

The proof is similar for Y in case (A) using the convexity at infinity of the drift, as we previously 
remarked. > 



Proposition 2.2 IfW satisfies (A) or W and V satisfy (A'), for all k E N there exists C{k) > 
such that for all 1 < i, j < N , 

supe(|X;'^ - ^P*^) < C{k) (1 + /xo ® /xo {\x - yp*^)) . 

Proof 

< We write the proof in the case V = (i.e (A)), the case (A') is similar. 

Recall that all particles are exchangeable. We may apply Ito's formula up to the exit time of a large 
ball (for the whole system) and then go to the limit in order to get, 

N N r-f N 



E ^ - = E J] - X^'^l^'^ + 2k{2k - 1) / E ^ - XJ" 

i,j=l i,j=l ^ i,j=l 

J (ywixi'"' - xi'"") - vwixi'"" - x]'^)) .{xi'^ - xl'^)|x^'^ - xi'^'^'^-'^ds. 



E 

N 



Denoting Ak{t) =EJ2^j^i l^t'^ - ^/'^P'', and using (2) we obtain 

Ak{t) < Nfio^ IJ.o(^\x-yf''^ +2k{2k-l + C) j Ak-i{s)ds - 2kX Ak{s)ds . 

Applying Gronwall's lemma and an easy induction we thus have that A^ (t) is finite for all t. Accord- 
ingly we may replace the pair of times {0,t) by {t,t + e) and prove that 1 1— Ak{t) is differentiable. 
Differentiating at time t yields 

A[{t) < -2\Ai{t) + 2N{d + C) . 

Gronwall's lemma yields the desired result for A; = 1. The proof follows, using this bound and an 
easy induction. > 



We are able now to generalize Proposition 2.1 and get uniform moment estimates of every order 
(under assumptions on the initial condition). 

Corollary 2.3 // W satisfies (A) or if W and V satisfy (A') (where m is defined), then for all 
t> Q and all k~^ 1 there exists a constant c{k) such that for all i 

supE\xi'^\^'' < c(fc)(l + ^o(k|^"'^)) . 



t>0 



IfW satisfies (A) a similar result holds with Y^'^ in place of Xl'^" . 
Proof 

<l Wc write the proof under {A'). Let Bk{s) = K{\Xl'^ \^). As in the previous propositions we shall 
use Ito's formula up to the stopping time T and then go to the limit. Using our hypotheses we get 
that for some nonnegative A' 

B2k{t) < B2k{0) + 2k{C + 2k- 1)E /* \Xi^^\^^-^ds - 2kX' E f \Xi'^f''ds 

Jo Jo 

i=i -^0 



,Ni2k-2 



ds. 



To bound the last term above, we use (A3) i.e. 

\vw{xi^^ - xi'^)\ < M(i + \xi'^ - xi'^D , 

and Holder inequality in order to obtain the following upper bound 



2(2fc-l)/2fc jy- 



E\X 



i,N\2k 



ri,N vi,N\2mk 



E(1 + |X]''^ -XJ- 



l/2fc 



ds . 



Now we may use exchangeability and Proposition 2.2 to obtain 

B2k{t) < B2k{0) + 2k{C + 2k-l)E r\Xi'^f''-'^ds-2kX'E TlXi'^^f^ds 

Jo Jo 

ft 2fc-l 

+c(fc)(i+/xo(ki''"'')'/''= y [E^ixi'^p*^)] ds. 

As usual, one can get that B2k is differentiable and satisfies 

B!,,it) < c{k)B2k-2{t) - 2kB2k{t) + c{k){l + ^io{\x?'^''))B2k{t)'-^ . 
Since, for every e > 0, it exists c such that jxp*^"^ < c + e|.Tp^ and a~2^ < 1 + a we thus obtain 

B'2k{t) < c{k){l + /Uo(|x|2™^)) - XB2k{t), 
for some A > and we can conclude using Gronwall one more time. > 



Remark 2.4 IfV is identically we may simply remark that 




and directly conclude with the help of Proposition 2.2, thanks to the convexity of x ^ |x|^. One 
could mimic the preceding proof to get the uniform estimate of the moments of the projected particle 
system under assumptions (A). 



2.2 Identification and existence of solutions of the nonlinear SDE 

We prove here how a tightness criterion may ensure that the empirical law of the particles converges 
to the solution of the nonlinear SDE (3) proving thus the existence of such a solution. 

Lemma 2.5 IfW satisfies (A) or W and V satisfy (A'), then for all s and t smaller than T, all 
k ^ 2 and all i 

E\xi'^ - X^'^l^'^ < C{k,T)il + Mo(kr^'"+^'=~^))) \t - s\i . 

Proof 

<3 Let < s < t < T and A; ^ 1. 



+ (2A;)E^* \\7V{Xii^)\\Xii^ - Xi'^f'^-^du 

Hence using (A3), Corollary 2.3 and Holder (with the best choice of exponents) we get that for 
k^l 

E\xi^^ - Xy^\^^ < C{k,T){l + fioi\xr^"'+^''~^^)) \t-s\. 

Plugging this estimate into the previous inequality and using Cauchy-Schwarz , we obtain the de- 
sired result. > 

It is well known that Lemma 2.5 implies that the sequence of the laws of (s —>■ xI'^)n defined on 
C([0, T], di?*^) is tight. In order to build a solution to the nonlinear SDE (3) wc may now follow 
some standard routine in mean-field particle systems. We here follow the one in [Mel95] Theorem 
4.1.4. Thanks to Proposition 4.2.2 in [Mel95] and to Lemma 2.5, the empirical measures ttat defined 
on r{V{C{[0,T],R'^))) by ttn = l/N ^l^i Sxi,N is tight too. According to the end of section 4.2 
in [Mel95], for any limit point vToo, any Q G P(C([0, T], M^)) IS TToQ ct.S. ct solution to the nonlinear 
martingale problem (up to time T). Actually on one hand the proof here is simpler since we have 
no jumps, but on the other hand the drifts are unbounded (but with polynomial growth) so that 
to justify passage to the limit we have to use the (uniform in N) moment estimates in Corollary 
2.3 and Proposition 2.2. Details are straightforward and left to the reader, who will check that a 
sufficient condition is //o(|a;p ) < +oo. 

But consider two solutions (X , Z ) of the nonlinear SDE (3), built with the same Brownian motion 
and the same initial condition, and introduce an independent copy (X' , Z'). Thanks to Proposition 



2.7 in the next subsection, we know that all processes have a finite second order moment. Hence, it 
holds 

A{t)=E[\Xt-Zt\^] = - I 2¥.[{yV{Xs)-VV{Zs)).{Xs-Zs)]ds 

Jo 

- [ 2E[{VW{Xs-X',)-VW{Zs-Zi)).{Xs-Zs)]ds. 
Jo 

But again since VW{—a) = — VVF(a) the later integral can be rewritten 

- /* E[{\7W{Xs - Xi) - \7W{Zs - Zi)). {{Xs - Xi) - (Z, - Zi))]ds . 
Jo 

Using (2) and the local Lipschitz property we see that -{VV{x) - W{y)).{x - y) < P\x - a 
similar result holding for W. Using (a + b)^ < 2a^ + 26^ and the fact that X — Z and X' — Z' have 
the same law we get 

A{t) <Sl3 A{s)ds, 
Jo 

so that A{t) = 0. Hence we have proved (strong) uniqueness for the nonlinear SDE. As for linear 
SDE, this notion of uniqueness implies uniqueness in law. 

As a byproduct, we obtain that tToo = 5q for some Q which is the unique solution of (3) and that 
TTN goes to 5q. 

What we have obtained is the following: there exists an unique probability measure Q defined on 
C([0,r],M'=') such that for all smooth / 



(a;, t) ^ f{ut) - f{uo) - f {Af{us) - VV ■ Vf{us) - VW * Q,(a;,) • V/(a;,)) ds 

Jo 

is a (5 martingale, with 

VW * Q,(a;,) = j VW{us - y)Qs{dy). 

Let H{x) = JVW{x-y)Qs{dy). Thanks to Corollary 2.3, H{x) < C(l + |x|"^) (provided Mo(|yr') 
is finite) and is local Lipschitz (thanks to (A3)). Thus the SDE 

dZt = V2dBt - VV{Zt) dt - H{Zt) dt 

has a strongly unique solution up to its explosion time. Since Q is a solution, it is the only one. The 

solution is thus non explosive, and Girsanov theory tells us that Q is absolutely continuous w.r.t. 
the Wiener measure, provided the drift is of finite energy (see e.g. [CL94] Proposition 2.3), i.e. 



J (|VF|2 + \H\^) (x) Qsidx) ds < +00. 



This condition is satisfied provided no (^\x\'^"^^^ < +oo. Thus for all s > 0, Qs is absolutely 
continuous w.r.t. Lebesgue measure, with a density Ug. 

Let us summarize our results. 

Theorem 2.6 Assume that W satisfies (A) or that W and V satisfy (A'). Assume in addition that 

^od^^l") < +00 for a = max(m(m + 3), 2m^). 

Then the nonlinear SDE (3) has an unique strong solution Q. 



Furthermore for all t > 0, Qt (the law of Q at time t) is absolutely continuous w.r.t. Lebesgue 
measure, with density ut satisfying for all T > 0, 

sup / \x\^''utix)dx < C(r)(l + //o(k|""'^)) , 

Q<t<Tj 

and t ^ Ut is a solution of (1). 

Finally t utdx is the unique solution of (1) among the set of continuous flows of measures t ^ Uf 
satisfying for all T > 0, 

J \x\^'^^dvtdt < +00 . 

The only thing it remains to prove is the last statement. Let t ^ be a solution of (1) satisfying 
the (finite energy) condition above. Then, according to Theorem 4.18 in [CL94], a solution of the 
(linear time inhomogeneous) SDE 

dZt = y/2dBt - VV{Zt)dt - {VW * ut){Zt)dt 

with initial law ^ug exists and furthermore its law at time t is given by Vf Hence Z is a solution of 
the nonlinear SDE. Uniqueness of the later implies that ut = utdx. 

About the literature. If existence and uniqueness have been extensively discussed in the frame- 
work of mean field interacting particle systems with bounded interactions, the case of unbounded 
interactions was not much studied. The above proof shows that convexity at infinity allows us to 
essentially mimic the bounded case, just changing the tightness criterion to be used, and using some 
recent aspects of stochastic calculus related to singular diffusion processes. 

If d = 1 another version of Theorem 2.6 is obtained in [BRTV98] Theorem 3.1 following a completely 
different way. One can notice that some moment condition on jiQ similar to ours is also required 
therein. Our condition (as well as the one in [BRTV98]) is certainly non sharp. If a large part of 
the method in [BRTV98] can be extended to the d dimensional case, some aspects (strongly using 
monotonicity) require additional work. 

For the nonlinear PDE, for d = 1 and m = 2 an existence and uniqueness result is stated in 
[BCCP98] p. 983. However the initial measure has to be absolutely continuous with a density 
satisfying J |x|^(i/xo < +oo (while we need a 10 instead of a 4), and uniqueness holds for classical 
(C^) solutions. 

The conclusion is that such results are certainly not useless. Since we shall need stronger moment 
assumptions in the sequel, we did not try to obtain the sharpest conditions in Theorem 2.6 (the 
interested reader may indeed remark that one need for the use of lemma 2.5 for tension only a 
dependence in time of the order \t — s\^^'' so that the condition on the initial measure may be 
weakened). 

2.3 Uniform moment control for the nonlinear SDE 

We start with the proof of the moment control for any solution of the nonlinear SDE. 

Proposition 2.7 Assume that W satisfies (A) or that W and V satisfy (A'). Then 

supE|Xt|2 </io(kP) + if. 

Proof 

< Let s ^ X'g SM independent copy s ^ Xg. Then 

E(|Xtp) = E(|Xop)-2fc/ ¥.[Xs-{VV{Xs) + VW *Us{Xs)))ds + 2dt 

Jo 

= E(|XoP)-2 / E{Xs-VV{Xs))ds- [ E{{Xs-Xl)-VW{Xs-Xl))ds + 2dt 
Jo Jo 



so that if we denote v{t) = E(|Xtp), we may first differentiate w.r.t. time as we did before, and 
tfien use the hypotheses. Remark that, since W is convex at infinity and symmetric, 

X ■ VW{x) > X ■ VW{0) + \w\xf -C^-C. 

It follows 

v'{t) < -2Xvv{t) + {2d + 3C) - 2\w {v{t) - E^{Xt)) ■ 

If (A') holds, \v > and Gronwall's lemma concludes the proof. If F = 0, \w > and we have 
assumed that E{Xt) = so that Gronwall's lemma concludes the proof. > 

Now, if we assume that existence and uniqueness hold for the non linear SDE, the results in the 
previous two subsections imply some control of the moment of the solution, just taking limits. But 
the previous a priori bound was necessary to complete the proof of Theorem 2.6. 



3 Propagation of chaos 

This section is devoted to the comparison of the behavior, for a fixed number of particles, of the 
difference between one particle and the solution of the non-linear SDE. Ideally, it will be uniform on 
time and will decrease quickly to as the number of particles increases. We will however see that to 
get such an estimate, we have to introduce a new convexity assumption, ensuring strict convexity 
except at some points. 

All results in this section are written under assumption (A) i.e. with a vanishing confinement 
potential V . Replacing Y by X, the same results hold under assumption (A'), modifying the proofs 
in the same way as we did for various statements in the previous section. 



3.1 A first control on the error 

We prove here a first uniform control on the mean square error between one particle and the 
solution of the nonlinear SDE. In the sequel we shall always make assumptions ensuring existence 
and uniqueness of strong solutions. Hence we may build solutions for (8) (or (7)) and (3) with the 
same Brownian motions and the same initial random variables (obtaining thus and i.i.d. sample 
{X%=i,...,N oiX). 

Theorem 3.1 Assume that W satisfies (A). Suppose that the law admits a large enough poly- 
nomial moment. Then there exists K > such that 

supEf|r/'^-Xif ) < i^. (9) 

t>o ^ ^ 

This result is of course not so good as it doesn't imply the so called uniform propagation of chaos, 
since the last bound does not tend to as A*" tends to infinity. However under the condition that 
W is the sum of an uniformly convex function and a Lipschitz compactly supported one, one may 
show the following bound, for some positive a 



E 



which gives us the non uniform propagation of chaos. The proof of such an estimate is the conjunc- 
tion of the proof below and standard estimates for bounded Lipschitz drifts for mean field particle 
systems (see again [Mel95]). 
Proof 
< Introduce 

1 ^ 



Then 



2 



= 2E(|y/'^-y/i2)+-E(|xip) 



since E(X|) = and the Xl 's arc independent. The second term is bounded by some K/N according 
to Proposition 2.7 so that we focus on the first term. But 



_ yhN 



so that 



with 



^ -^0 j=l 

^ _ 9 ^ ft 



^,,■(5) = {vw{Y:'''-Yi''')-vw{Y:'''-Yj'''))-{Y:'''-Y:'''), 

B,,{s) = {VW{Xl-Xi)-VW*Us{Xi)).{Yt'' -Yt""), 

Cij{s) = vw*u,(xi)).(y;'^-y;.^). 



Let us deal with the first term, 

N , N 



TV 



using (2) and properties of vectors in Ai. 
Remark now that 

N ( N 

-Ej^Bijis) < (E|y;'^-y;'^|2)^/ME|E(VM^(^i-^i)-vw*^z,(xi)'2 
< cv^(E|y;'^-y;'^p)'/' 

by Cauchy-Schwarz inequahty, the polynomial growth of VW, the controls of moments previously 
established and the key remark 

E{VW{Xi - Xi) - VW * Us{Xi)) = 0. 

The last term involving Cij{s) is controlled using the same tools. We finally get, defining 

a{s)=E(\Y:^'^ -Y:'^^) 



and reasoning as we did in the previous section, the differential inequahty 



a'{s) < -2Xa{s) + -j=^/a(s) + C/2. 

Now we may use y/a{s)/^/N < l/2(e~^-^ + ea{s)) and choose £ = 2A/c. The previous inequality 
becomes a classical 

a'{s) < -Xa{s) + (c/eiV) + C/2 
for which we can use Gronwall's lemma. > 

Remark. When C = we get the uniform (in time) propagation of chaos result obtained in 
[Mal03] with rate 1/N. Of course in our case (convexity at infinity only) such an uniform propagation 
of chaos cannot be expected (more precisely cannot be shown with usual tools). Indeed, the collision 
potential W obliges the particle to mainly stay in the same region, where there is no more attraction 
between particles. 

3.2 Uniform propagation of chaos 

In view of Theorem 3.1 and its proof, we have to reinforce the convexity assumption in order to 
prove the uniform propagation of chaos phenomenon. 

The condition (6) in the introduction and inspired by the work of Carrillo-McCann-Villani [CMCV06] 
is recalled below: 

we say that condition C(A, a) holds if there exist A,a> such that for any < e < 1, 

Vx, y e R'^, {x - y).{VW{x) - VW{y)) > ^e"(|a; - y\^ - e^). 
We can now prove 

Theorem 3.2 Assume that W satisfies C(A, a) and (A). Suppose that the law fj,o has a large 
enough polynomial moment. Then there exists K > such that 

supE(|y/'^-X,^p)<^. (10) 

Proof 

< The proof follows the same lines than the one of Theorem 3.1. First recall that 



E(|y/'^-x,^|2) < 2E(|y/'^-y/|2) + |E(|x,| 



and the second term is of a better order thanks to Proposition 2.7 so that we focus on the first term. 
For this term the only modification is the control of Aij(s) where we replace convexity at infinity 
by (6). This yields 

N / ^ \ 

i,j=l \ i=l J 

SO that the differential inequality satisfied by a{t) becomes 



a\s) < -2^e"(a(s) - e^) + -^^/^. 



Applying Theorem 3.1 wc know that a{s) < K for some K > 1. 
We may take e = \/a{s)/2'/K < 1 and get 



a' is) < -Jaisf+''''^ + -^.Ms) 



/3'(s) + (J/2)/3i+"(s)< 



so that 

Pis) < C/Af^/2(l+a) 

for any s such that /3'(s) > 0. Since /3(0) = it easily fohows that (3{s) < C /N^/'^^^'^"^ everywhere, 
hence the result. I> 



Remark. If (A') holds we have to assume that V satisfies C(A, a) and W is convex. 



4 Ergodicity of the nonlinear SDE 

4.1 Drift condition and existence of a stationary measure 

In fact, as long as we consider only the particle system, the condition (2) is sufficient to ensure the 
ergodicity, and even exponential ergodicity, of the particle system. Indeed we have that (^i'^)i<j<Ar 
is strongly Feller aperiodic and we can use Down-Meyn-Tweedie's drift condition [DMT95]. This 
condition is the following: there exists some $ > 1 with compact level sets such that for some 
positive A, h 

< -A$ + b 

where C is the generator of the particle system. It is trivially verified here with the function 
<I>(x) = + 1 according to (2). It of course implies the existence of an invariant measure tt^ 
and that there exists 6{N) > such that denoting P/^ the semigroup associated with the particle 
system 

||P,^(a;,.)-vr^(-)llrv<ce-*^W<&(x). 

As a byproduct, we get the existence of an invariant measure for the nonlinear SDE (hence an 
equilibrium for the nonlinear PDE), provided we have uniform (in time) propagation of chaos. 
Indeed the previous inequality together with uniform propagation of chaos show that the family 
{ut)t>o is a Cauchy family in L^(da;), which is a complete space. 

Such a result is however not useful to control the convergence to equilibrium for the nonlinear SDE 
(or even the nonlinear PDE). Indeed there is no close form for S(N) which heavily depends on A*", so 
that even the uniform propagation of chaos cannot give good estimations for the rate of convergence 
of (Xt) to its invariant measure. Note that even a strict convexity condition does not seem to be 
useful for such an approach. 



4.2 Logarithmic Sobolev inequality 

Under a global uniform convexity condition (i.e. C = in (2)), Malrieu proves in [Mal03] that the 
invariant measure of the centered particles system i.e. 

/ 1 ^ \ 

where dy is Lebesgue measure on the hyperplane Ai, satisfies a logarithmic Sobolev inequality with 
a constant 1/A independent of the dimension using the famous Bakry-Emery criterion. So that if 
uf^ stands for the law of (l^*'^)i<i<Ar we have that 



Projecting on the first coordinate, using a T2 inequality and, on one hand the dimension free 
exponential convergence, and on the other hand the uniform (in time) propagation of chaos, he 
essentially proves the following 

W2(ut,Hoo) < Ce-^VEnt(no| 

This result was obtained by Carrillo-McCann-Villani [CMCV03], and improved in Carrillo-McCann- 

Villani [CMCV06] where the authors replace the square root of the initial entropy 1)}- Tr2(uO)^too) 

(recall that a log-Sobolev inequality implies a T2 inequality, i.e. ^2(^*05^^00) ^ CiyEnt('Uo|'"oo))- 

It is of course tempting to use the same approach in the non strictly convex case. 

The first point is that under assumption C(A,a), we may no longer use Bakry-Emery condition 

for the measure as the Hessian may degenerate. We are thus obliged to use some perturbation 

argument. 

To fix ideas and to simplify the presentation, let us restrict ourselves to a polynomial potential, 
namely W(x) = Ixp"*"". Consider first a crude attempt: let /3(x) = aa;^/(l + x"^) with some positive 
a: (3 is bounded, and has positive curvature 2a at 0. Choose now a = 1/N and consider the 
perturbed measure 

vg{dy) = Z^^exp |^/3(|y, - y,|)j u^{dy). 

By arguments developed in Malrieu [Mal03] one may then easily establish that has positive 
curvature bounded below by Thus it satisfies a logarithmic Sobolev inequality with constant 

eA'^ remarking that the perturbation potential has oscillations bounded by one and applying the 
standard Holley-Stroock perturbation argument. We may then use the entropic convergence asso- 
ciated to such a logarithmic Sobolev inequality, the T2 inequality projected on the first coordinate 
and uniform propagation of chaos (Theorem 3.2) to get 

/ \ 1/2 

W2{ut, Uoo) < + e-*/^"^V2eA^Ent(uo|?Xoo)- 

Optimizing in ensures that 

K' \ogm t 

W2{ut,Uoo) < -n— 

t2(l+a) 

which is quite a bad estimate of the speed since Carrillo-McCann-Villani obtained in such EL CclSG cl 
speed i-(V«) ( 

we shall recover this speed in the next subsection). 
Note finally that this approach via logarithmic Sobolev inequality has two major drawbacks: on one 
hand the right hand side depends on the entropy of the initial data, on the other hand it cannot be 
used to compare the behavior of two solutions starting from different initial measures. We will see 
in the next section how to solve both these problems by a direct simple approach. 

4.3 A direct control of the L^-Wasserstein distance 

The goal of this section is to show that even if the logarithmic Sobolev inequality cannot give useful 
result, one may use a direct approach to obtain a (uniform in N) control of the Wasserstein distance 
between two solutions of the particle system starting from different points. So, combining this with 
the uniform in time propagation of chaos, we recover in an elegant and simple way the results of 
Carrillo-McCann-Villani [CMCV06] . 

Theorem 4.1 Assume that W satisfies (A) and the convexity condition C{A,a). Let Uf and vt be 

the unique solutions of the nonlinear PDE with initial conditions respectively uq and vq. We assume 
for simplicity that both uq and vq have an exponential moment (or a large enough polynomial moment 
in order to ensure existence and uniqueness). 



Then t i— W2{ut,vt) is non-increasing. Fkirthermore there exists ti < (2^+"/3) log(W|(uo, fo))/^ 
such that 

Wi{ut,vt) < e-(3^/2'+")*VF|K,t;o) ^/^<^l, (11) 
< + ^ (a/(2 + a))^+"/2(i_t,)^"'/" ift>ti. 

In particular, ift > (1 + ?7)(2^+"/3) \og{W2{uo-,Vo))/A for some rj > 0, one has 

W^{ut,vt)<(l + iA/{l + r,)) (a/(2 + a))^+"/^ t)''^" . 
In addition for all t >Q one has 

Under an uniform convexity condition for W (i.e. C{A,Q)), the convergence is exponential i.e. 

W2{uuvt) < Ce-'^'W2{uo,vo). 

Note that the final assertion completes the result of Malrieu [MaI03] who cannot compare via 
logarithmic Sobolev inequalities the distance between two solutions. 

Remark also that the bound obtained for the control in Wasserstein distance in the convergence 
towards the stationary measure is also better as it is now ^2(^*05 Uoo) which controls the decay and 
which is smaller (since a T2 inequality holds) than C ^yEnt{uo\u^) . 
This will be explained later. 

As will be seen from the proof, there is no use here of heavy (but sharp) technology as optimal 
transport in length space as developed by Carrillo-McCann-Villani [CMCV06] and thus we believe 
that such an approach will give good results for other models (reinforced diffusion,...). 
Proof 

<1 In what follows we shall always consider coupling consisting in picking the same Brownian motions 
(B^) for the two particle systems, and only choose an ad-hoc coupling for the initial random variables. 
Denote ul'^ the law of the first particle of the centered particle system (Y) starting with initial law 
uq and consider v]' the law of the first particle of the centered particle system starting with initial 
law vq. Let us first remark that for each coupling go of uq and vq 

where (F*'^) is given by the centered particle system where each particle starts with law uq , (F'*'^) 
with measure vq and the initial law of the pairs {Y^^, Y^'^) are given by independent copies of ^o- 
The subscript for the expectation is related to the initial law. 
We then get by Ito's formula and symmetry (starting now from points y and y') 

i=i i,j •'^ 



with 

We may now differentiate in time and tfien use condition C{A, a) to get 

\i=l / i,j 



-Ae'^^y,y'T.(\yt'"'-y"'"'\"-^") 



which gives us by exchangeabihty 

d 



dt-y^yiiy:'"' - y^'^'i') < -Ae-n\yr - yl^'^'f - e'). 



Denote ^{t) = Ey^y> - 1"/''^^) , we then obtain for all e < 1 

This inequality for e going to 0, implies ^'{t) < 0, i.e. ^ is non- increasing. As a byproduct we get 
that t W2{ul''^ ,vl'^) is also non-increasing. Indeed if we choose go as the optimal coupling for 
the quadratic cost, it holds 

and we get the result just shifting the initial time. 

Now we can separate two cases: either ^(i) > 1 or ^(i) < 1. Note that there exists some ti > such 
that the first case holds for t < ti and the second one for t>t\. 

li t < ti we may choose any e, for instance here we choose e = 1/2 and obtain (since ^{t) > 1, 
£2 = 1/4 < C(t)/4) 

^'(t) < -A{a)m 
with A{a) = (3A/4) (1/2)", which gives by Gronwall's lemma 

e(t)<e-^(")*e(o). 

Choosing again the optimal coupling qq we obtain the first part of the result and the fact that 
ti<log{Wi{uo,vo))/A{a). 

For the second part, if t > ti choose = a^{t)/{a + 2) to get 

a/2 



Integrating this differential inequality we get 

e(^)<(l + s(a)(^-^l)^^/^ 

with B{a) = A (a/(2 + a))^+°/l Our choice of e gives the optimal constant (for this method). 

Of course writing C{A,a) as we did is a little bit artificial, and if we want a more homogeneous 
estimate we may remark that 0(^4,0;) implies 

(x - y) • {VW{x) - VW{y)) > ^(6/1^2(^0, vo)ri\x - - e^) 



for e < W2{uo,vo). So we may always choose = a^{t)/{a + 2) and get for all t, 



We use then the uniform propagation of chaos property to transfer the inequality from the particles 
to the solutions of the nonlinear SDE: 

W2{ut, Vt) < W2{ut, ul'"") + W2{ul''', vl'"") + W2{vt, vl^""), 

and take the previous (uniform in N) estimation for the middle term and the uniform in time esti- 
mation for the first and third term and let N go to infinity. Of course if ti depends on N, its bound 
does not, so that we may find a converging subsequence and get the result. > 

Remarks. We have seen in subsection 4.1 that there exists an invariant measure Uoo for the 
nonlinear PDE. The bound for W2{ut,vt) obtained by introducing the particle system of order N, 
and then choosing choosing vq = allows us to prove this existence too (using completeness of the 
Wasserstein distance). In addition the decay of the Wasscrstein distance trivially implies uniqueness 
of tioo) at least in the set of measures having some large enough polynomial moment. 

Hence, we recover in the previous result the asymptotic rate of convergence to equilibrium obtained 
by Carrillo-McCann-Villani. But it seems that the result improves upon theirs for small times where 
we are able to describe some exponential decay. However the value of ti is not explicit, so that for 
practical issues, this initial exponential decay is not really tractable. 

We choose here to develop the point of view of the particles; prove the decay of the Wasserstein 
distance for the particles system and then transfer this decay to the solution of the nonlinear SDE 
via uniform propagation of chaos. It is however easy to develop the same line of proof directly for 
the control of the Wasserstein distance of two solutions of the nonlinear SDE. 

Note finally that we did not consider here coefficient diffusion other than constant, as granular 
media equation is formulated with a constant one. It would however not be difficult to introduce a 
condition enabling to obtain the same decay. Indeed consider 

dZt = b{Zt)dt + a{Zt)dWt 

where Wt is the usual Wiener process in W^, and assume that for all < e < 1 there exists positive 
A and a such that 

{x-y)- {h{x) - h{y)) + \tT{{a{x) - a{y)){a{x) - a{y)f) < -Ae^{\\x - yf - e^) 

then we easily derive the same polynomial decay for the Wasserstein distance of two solutions Zf and 
Zt with different initial conditions. Note that this type of condition hold with possible degeneracy 
of the diffusion coefficient. 

Again in this section we should assume that (A') holds, V satisfies C(A, a) and W is convex. 

5 Concentration inequality 

The main goal of this section is to complete the results on convergence of the particle system 
and of the nonlinear system by providing a (non asymptotic) deviation inequality. This inequality 
will be written for additive functionals of the particles and then allows us to estimate integrals 
with respect to the stationary measure of the nonlinear PDE. Once gain it is interesting to get 
uniform in times estimation to be able to simulate at fixed time the particles and use them for the 
evaluation with no loss at each time of the constant in the concentration. Therefore we first prove 
in the general framework an uniform Ti inequality for solution of SDE subject to some convexity 
at infinity condition. We then show how to use them on our example. 



5.1 Uniform transportation cost inequality under convex at infinity condition 

Let X be the solution of the following stochastic differential equation: 

dXt = dBt + b{Xt)dt, (12) 

where i? is a standard Brownian motion on and b a smooth function from to W^. 
In the case when b{x) = x, X is the well-known Ornstein-Uhlenbeck process and 



sup IE 



J\X^-Y,-f 



< +00, 



(where X^ and are two independent copies starting both at x) is finite if and and only ii S < 1/2, 
since the law of Xf is the Gaussian measure 

Caw{X^) = M xe~\ 



Our aim is to extend this assertions to the case when the drift is confining only outside a compact 
set, and for non-constant diffusion coefficient, so that X is the solution of 

dXt = a{t)dBt + b{Xt)dt. (13) 

We use here the formalism of transportation cost inequalities in W\ distance for which a practical 
criterion based on the integrability of the exponential of the square of the distance is sufficient 
[DGW04] (and Bolley-Villani [BV05] or Gozlan [Goz05] for a better evaluation of the constant) and 
which implies interesting deviation inequalities. 



Proposition 5.1 Suppose that there exist A, A > and C such that, for every x,y & 

{x-y)- {b{x) ■ 
then for any S < X/2A and any x G 



{x-y)-{b{x)-b{y))<-X\x-yf + C, Mhs < A (14) 



supE 



< 1 + (Ad + C + i)eS(^d+c+mx-2SA) ^ 



where X^ and are two independent copies of (12) starting at x. 

So the law of Xf satisfies a Ti inequality with a constant £ independent of time and initial position. 
Recall that a Ti inequality with constant £ for a measure n reads as: 



Wiiiy,!!) < V^Ent(z^ 

for any v. 
Proof 

< The process Z = X — Y is a solution of 

dZt = cj{Xt)dWt' - a{Yt)dW^ + ib{Xt) - b{Yt))dt, 

with initial condition Zq = 0, where and are two independent Brownian motions. Ito's 
formula ensures that 

^s\x,-Y,f _ ^s\x.-u^ < 26l\x^-Y^).{b{X^)-b{Y^))e'\'-\'du 

ft 2 
+2dA J (d + 2(5|X„ -y„p)e^l^"l du + Mt-Ms 



< 26 J {Ad + C+ {2AS - X)\Xu - y„|^)e'^l^"l'du + iW* - M^, 



where M is a local martingale with quadratic variation given by 

{M)^ = 8.4(52 f\X^ - Y^\^ e^\^--^^\^ du. 
Jo 

For every S < X/2A, let us denote by K the quantity 



K 



Ad + C + 1 
X-2SA ■ 



For every x G M , 

{Ad + C-{X- 2M)|a;|2)e^l^l' < (Ad + C - (A - 2M)|i^|2)e^l^l'l(|^|^K) + {Ad + C)e^^^ 

< -e''l^l' + (Ad + C+l)e'^-^'. 

For R> 0, introduce Tr the first time when Z exits the ball of radius R and define an the function 
defined by 

Then 

adt) - 1 < 25(Ad + C + l)e^^\, 
so that we may let R go to infinity, and have shown that the exponential moment is finite. 



Defining now a{t) = E 



=-51^*1' 



, one gets, for < s < i. 

As a consequence, a satisfies the following differential inequality 

a'(i) < -25a{t) + 25{Ad + C + l)e^^\ 
Denoting by (3{t) = a{t) exp(25t), this implies that, 

I3'{t) < 25{Ad + C + l)e^^' e'^^\ 

and then. 



u) du. 



1 . 



Pit) < P{0) + {Ad + C + l)e^^' (e^^* 

Besides, X and Y have the same initial condition so /3(0) = a(0) = 1. As a conclusion, we get for 
all t ^ 0, 

a{t) < + {Ad + C + l)e^^\ 

which achieves the first part of the proof. The final statement follows from [DGW04] Theorem 3.1. 
An explicit expression of C is derived in [Goz05] chapter VII or in BoUey-Villani [BV05]. > 

An important consequence of a transportation inequality is that we easily obtain deviation inequality 
for Lipschitz functions. Indeed for all Lipschitz functions F with ||-F||Lip < 1 and all positive r 



Remark that the previous proof extends to the case when the initial law /xq satisfies 

e^\''-y\\Q{dx) noidy) < +oo 
(that is 1^0 satisfies a Ti inequality) just choosing two independent variables Xq and Yq of law //q- 



/ 



5.2 A concentration inequality for the stationary measure 



We may now recombine results inherited from the previous sections to get an useful inequality for 
the evaluation of / fduoo when / is a Lipschitz function. Prom the previous section, if (A) holds 
and /io satisfies a Ti inequality, the particle system satisfies a Ti inequality with constant €iV (for 
some <t independent of time) , which thus leads to the following: 
for / Lipschitz (in M"*) with ||/||Lip < 1, for all positive r and all t 



\ k=l 



>r\< e-^'-'/^. 



Remark now that for a Lipschitz function satisfying ||/||Lip < 1, if the convexity condition C(A, a) 
holds, the uniform propagation of chaos of Theorem 3.2 tells us 



|E/(X,^'^) - J f{y)ut{y)dy\ < wM'"" ,^t) < W^iu]'"' ,ut) < 
so that for all r > ( — ^ 

^ hv E /(^'''^) - / fiy>t{y)dy >r-[^ 



K 



1 



1/2 



1/2 



We may now use convergence in W2 distance of the solutions of the nonlinear SDE towards the 
stationary measure given in Theorem 4.1 to get 

/ \ 1/2 

Proposition 5.2 For all r > [ — ^ + JM), 

^ E /(^'''') - / fiy>Mdy > r - (-|-) - < e-^'^Vc^, (15) 

where (3{t) is one of the functions governing the decay of W2{ut,Uoo) described in Theorem 4-1- 

It is quite hard to imagine to extend to our case uniform result over Lipschitz function (i.e. deviation 

of the Wi distance between the empirical law of the particles and the stationary measure) as in 
BoUey-Guillin-Villani [BGV05] as it requires a dynamic coupling which can be achieved only for 
potential W whose Hessian is bounded (in the sense of matrix). 

Let us finally note than one can use an Euler-Maryama scheme preserving square exponential inte- 
grability and with good stability property [LMS06, Thl and Th 4] to simulate the particle system 
leading to the same concentration inequality. Note that the recurrence property needed for the 
stability of this adaptive scheme in [LMS06] is exactly our condition of convexity at infinity. 
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